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Dynamical magnetic and nuclear polarization in complex spin systems is discussed on the example 
of transfer of spin from exciton to the central spin of magnetic impurity in a quantum dot in the 
presence of a finite number of nuclear spins. The exciton is described in terms of the electron 
and heavy hole spins interacting via exchange interaction with magnetic impurity, via hypeprfine 
interaction with a finite number of nuclear spins and via dipole interaction with photons. The 
time-evolution of the exciton, magnetic impurity and nuclear spins is calculated exactly between 
quantum jumps corresponding to exciton radiative recombination. The collapse of the wavefunction 
and the refilling of the quantum dot with new spin polarized exciton is shown to lead to build up 
of magnetization of the magnetic impurity as well as nuclear spin polarization. The competition 
between electron spin transfer to magnetic impurity and to nuclear spins simultaneous with the 
creation of dark excitons is elucidated. The technique presented here opens up the possibility of 
studying optically induced Dynamical Magnetic and Nuclear Polarization in Complex Spin Systems. 

PACS numbers: 75.50.-y,75.50.Pp,85.75.-d 



I. INTRODUCTION 

There is currently interest in developing means of local- 
izing and controlling complex spin systems in soli d st ated 
This includes electron and/or hole spins in gatecP^, self- 
assemblecP, nanocrystaP and carbon nanotube quantum 
dotsP, nitrogen vacancies in diamoncP and magnetic im- 
purities in II- 

v pm and in .^SESl_ The complex spin 
systems involved include heavy valence holes with spin 
J = 3/2, nitrogen vacancies with spin M = 1, half- 
filled shell electrons of mangan Mn 2+ impurity atom 
with spin M = 5/2 or Mn 3+ atom with M = 3/2 in II- 
VI semiconductor quantum dots, or strongly coupled va- 
lence hole-Mn atom in InAs/GaAs quantum dots. Exten- 
sive theoretical studies have been carried out, predicting 
rich phase diagram for these system y^H For NV cen- 
ters in diamond, carbon nanotube based quantum dots 
and magnetic impurities in II- VI semiconductor quantum 
dots the complex spin system interacts with only a finite 
number of nuclear spins. The controlling of magnetiza- 
tion of complex spin systems is often carried out optically 
and involves transfer of photon angular momentum into 
exciton spin, a nd exc iton spin into the spin of the com- 
plex spin systerrP^ESl. This dynamical magnetic polariza- 
tion (DMP) process is decohered by photon and nuclear 
spin baths. Recently, first optical experiments on single 
magnetic impurities in II- VI quantum dots measured the 
dynamic evolution of the magnetization proces jHEH l w jth 
theoretical models of DMP based on rate equation d 24 * 25 l 

In this work, we develop a microscopic theory of opti- 
cally driven dynamical magnetic polarization of complex 
spin systems. The theory describes the transfer of spin 
from exciton to the central spin of magnetic impurity 
in a quantum dot in the presence of a finite number of 



nuclear spins using quantum jump approaclP3ES ^p^g 
exciton is described in terms of the electron and heavy 
hole spins interacting via exchange interaction with mag- 
netic impurity, via hypeprfine interaction with a finite 
number of nuclear spins and via dipole interaction with 
photons. The time-evolution of the exciton, magnetic 
impurity and nuclear spins is calculated exactly between 
quantum jumps corresponding to exciton radiative re- 
combination. The collapse of the wavefunction and the 
refilling of the quantum dot with new spin polarized exci- 
ton as in recent experiment by Goryca et alP is shown to 
lead to build up of magnetization of the magnetic impu- 
rity as well as nuclear spin polarization. The competition 
between electron spin transfer to magnetic impurity and 
to nuclear spins simultaneous with the creation of dark 
excitons is elucidated. 

The paper is organized as follows. In Section In] we 



describe our model. Section III describes quantum jump 
approach to time evolution of our system for a system of 
single MI, single exciton with no nuclear spin. Section |TV| 
contains quantum jump approach and the dynamical evo- 
lution of MI by train of excitons in the presence of nuclear 
spins. In sections[V|and [Vl] we presents numerical results, 
discussion, conclusion and the summary. 



II. THE MODEL 

We consider a semiconductor QD containing a complex 
spin system M, e.g., magnetic ion (MI), coupled with 
few nuclear spins of the host material, as in, e.g., CdTe 
quantum dots. The quantum dot with MI is attached to a 
smaller quantum dot with no MI where the electrons and 
valence holes with definite spin polarization are generated 
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FIG. 1: Schematic representation of the DMP process. 



optically by circularly polarized light. This is illustrated 
in Fig. la where circles describe quantum dots, blue arrow 
corresponds to electron spin S z — +1/2 and white arrow 
corresponds to heavy hole spin J z = —3/2 in the small 
dot. The large dot contains a randomly oriented complex 
spin M, represented by a magenta arrow, and a number 
of randomly oriented nuclear spins represented by small 
arrows. The DMP process starts with transfer of spin 
polarised exciton from a small QD to a large QD, Fig. lb. 
As a result of interactions resulting in the flipping of elec- 
tron, MI and nuclear spins the wavefunction of the large 
dot evolves in time and becomes a linear combination of 
bright and dark exciton states, Fig.lc. During this pro- 
cess the small dot is refilled with spin polarized exciton. 
Simultaneously in time the bright exciton decays due to 
interaction with the photon field, and at some point in 
time, it recombines, photon is emitted, and the quantum 
jump takes place. The states of the magnetic ion and 
nuclear spins are modified, polarization is increased and 
the large dot is refilled with spin polarized exciton and 
the DMP process continues. 

We now quantitatively describe the DMP process. We 
start with the Hamiltonian describing the quantum dot 
coupled with the photon bath H = Hqd + -f/phQD + 



i/ph . Here H p h is the photon Hamiltonian, HphQD is 
the Hamiltonian describing coupling of photons with the 
exciton in a QD and Hqjj is the QD hamiltonian. 

The QD Hamiltonian describes exciton X coupled with 
the magnetic moment of the complex spin system MI and 
with nuclear spins / given by Hqb = -£f m + H x + H XI[L + 
H n + H xn + H mn . Here H x describes the exciton, H m 
describes the complex spin system MI and the remain- 
ing terms in Hqo represent X-MI, X-I and MI-I nuclear 
spins exchange couplings. The exciton Hamiltonian de- 
scribes the low energy quadruplet |5*, J) characterized by 
quantum spin numbers of an electron, S ~ ±1/2, and 
a heavy hole, J = ±3/2 in the QD. The complex spin 
system is described by a total spin M = J2i=i ^i, where 
N is the number of spins it = 1/2 building up the MI sys- 
tem, and H m = J2i<j Jij^i - Uj +DM% with Jij exchange 
matrix elements building the total spin M. In quantum 
dots one often includes strain field D leading to splitting 
of the different M z levels^. Similarly I = Yli^i 1% where 
Ni, is total number of nuclear spins. 

We assume that exchange coupling constants of MI 
spins with the environment (including excitons and 
nuclear-spins) are identical. Hence the full QD Hamil- 
tonian can be written as^ 

H QD = H m + H X + J hm J z M z ~J cm S-M + AI- M 

N b 

+ ^[Jnc-fn ' & + JnkIz,nJz] + ^ J ™ ,7 » ' ^' ' W 

n n,n' 

The exciton Hamiltonian iJ x = AqS z J z + Ai(S f+ J _ + 
S~ J + ) describes splitting Ao between the low energy 
dark exciton doublet | t,1r) = I + 1/2, +3/2), | |,4) = 
| — 1/2,-3/2) with total angular momentum j z = ±2 
along quantization-axis, z, and higher energy bright ex- 
citon doublet | = | - 1/2, +3/2), | t,JJ-) = | - 
3/2, +1/2) with j z = ±1. Here t / I and ft / 
represent spin of electron and hole^. The bright exci- 
ton doublet is split by the anisotropic electron-hole ex- 
change interaction characterized by parameter Ai which 
measures the splitting of the two bright exciton states 
| + 1/2,-3/2), | - 1/2, +3/2). Ai is zero for cylindri- 
cal quantum dots and the two bright exciton states cor- 
respond to circular photon polarization. The exciton- 
MI coupling in Eq.l is given as a sum of the ferromag- 
netic Heisenberg electron-MI exchange H em — ~J cm S-M 
and anti-ferromagnetic Ising exchange interaction -ffhM = 
+J\imJzM z ^ Only electron-MI interaction is responsi- 
ble for the e-MI spin flip-flop process. The interaction 
of complex spin MI with nuclear spin associated with 
the spin complex is denoted here by Hmi = AIm ■ M. 
This interaction might, for example, describe coupling 
of manganese d-shell electron spins with manganese ion 
nuclear spirpSl. With hole spin strongly aligned along 
the growth z direction the coupling of electron and hole 
spins to surrounding nuclear spins of isotopes of the 
QD and barrier material with finite nuclear spin reads 

Enl^ ' S + J n lJ z ,nJz} + E„,n' Jim' In ' In' where N b 
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FIG. 2: Time evolution of density matrix with the initial 
condition pwit = 0) = p22(i = 0) = 0.5 and Pij{t = 0) = 
for other i and j's, corresponding to initial random state of 
MI, are shown in (a), (c), and (d) for 8 — 0,1,5 meV. In 
(b) the expectation value of MI spin for 5 = is shown. 
Here |1> = \X b ,t), |2> = \X b ,l), |3> = \X d ,t), |4> = |X d ,|>, 
5} = |0,t), 1 6) = |0,4). The population of vacuum can be 
calculated by p V acuum = P55 + pm- The elements of density 
matrix, not plotted in this figure, are all identical to zero. 



is the number of nuclear spins in the QD and the last 
term describes nuclear spin interaction. We note that for 
isotropic QD the long range e-h exchange Ai is zero and 
the heavy hole spin J z = ±3/2 is preserved. 



III. SINGLE MI, SINGLE EXCITON AND NO 
NUCLEAR SPIN 

We start our discussion of DMP by discussing time 
evolution of magnetization of MI interacting with an ex- 
citon in the absence of nuclear spins. To focus on quan- 
tum dynamics we consider the simplest complex spin sys- 
tem, M = 1/2, with just two states | t) = \M Z = 1/2) 
and I I) = \M Z = —1/2) and Hamiltonian Hqd — 
H x + H xm where H Km = -J em S e ■ M + JhmS z ,hM z . In 
the basis of a single exciton-MI complex |1) = |X&,t), 
|2) = \X b ,i), |3) = \X d ,t), |4) = \X d ,i), |5) = |0,t 
), |6) = |0, 1), the exciton Hamiltonian is diagonal 
H x = diag(E b ,E b ,E d ,E d ,0, 0) where E b , and E d are 
the energy of bright and dark excitons measured rela- 
tive to the vacuum. In the basis of {|1), |4)}, {|2), |3)}, 
and {|5),|6)}, the X-MI Hamiltonian is block-diagonal 
H X m — flxro.l © H xm ,2 © H xm , 3 where 7J xm ,l = (—Jem + 

(Jem Jhm)/4 Jem/2 

Jem/ 2 (Jem Jhm)/4 

and i? xm ,3 = respectively. Here 1 is a 2 x 2 unit ma- 
trix. 



J hm )/4 1, H x 



The off-diagonal elements of H xm _2 describe mix- 
ing of X b and X d via spin-1/2 MI, hence \ip(t)) = 
C bti {t) exp(-iE b t/h)\X b , |) + C rfit (t) cxp(~iE d t/h)\X d , f 
) with initial condition \ip(t = 0)) = |Xb,4-) is one of 
the solutions of the time-dependent Schrbdinger equa- 
tion —ih-g^\ip(t)) — HQu\4>{t)) . This state describes a 
coherent Rabi-oscillations between bright and dark exci- 
tons due to MI spin flip-flop. 

Note that in \ip(t)) there is no mixing with the vac- 
uum, |0, M z = ±1/2), unless we take into account the 
coupling of bright-exciton with radiation field. In the in- 
teraction and rotating wave approximation the electron- 
photon coupling is described by the Hamiltonian that 
does not directly change the state of MI 



H phQ u(t) = KY,9%[b\%M z ){X b ,M z 

k,M z 

+bj;\X b , M z }(0, M z \ e +l ^ b - UJ ^ t 



i(u b —uiT-)t 



(2) 



where bt and 6g are creation and annihilation operators 
of photon with one specific circular polarization, g^, and 
are the photon-X coupling constant and photon fre- 
quency, and ui b = E b /h. The equation of motion of the 
QD density matrix, p, coupled with thermal bath of pho- 
tons can be calculated after tracing over photons degrees 
of freedom. Assuming that photons are in thermal equi- 
librium and are weakly coupled with excitons in QDs, 
the equation of motion for exciton density matrix, p, can 
be calculated perturbatively. Up to the second order of 
perturbation, it is straightforward to show thatP^l 

'''' = - l -[H QD ,p]-^-n B Y i {%M z ){Q,M z \ P 



dt 



(3) 



-2\X b , M z ) (0, M z \p\0, M z ) (X b , M z \ 
+p\0,M z ){0,M z \) 

-\(n B + l)J2(\Xb,M z )(X b ,M z \p 

— 2]0, M z ) (X b , M z \p\X b , M z ) (0, M z \ 
+p\X b ,M z ){X b ,M z \), 

where n B — l/(e hu / kBT — 1) is Bose-Einstein distribution 

function, and T = g^J, is the transition rate for the 
spontaneous emission of photons. £ is the dipole moment 
matrix element. Note that in Eq.([3]), vacuum can be 
considered as a shelving-state. 

The numerical solutions of Eq. ([3| at zero-temperature 
(n B = 0) are shown in Fig. [5] for a QD with E d — 2 eV 
and 5 = E b — E d — 0, 1, 5 meV. Here we used J om = 1 
meV and Jh m = 4 meV. The initial state of MI is com- 
pletely uncorrelated with half of the spins populated in 
up-direction. As it is shown, because of the coupling with 
the bath of photons, bright-exciton decays to vacuum 
without flipping the MI spin and mixing with X d , e.g., 
\X b ,M z ) -> \0,M Z ). In Fig.[2j we find that p n = e" rt /2 



and P55 = (1 



-n 



)/2 fit perfectly the numerical so- 



lution of pn and P55 for all Ss. The decay channel of 
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FIG. 3: Time evolution of S z of a train of injected photo- 
electrons inside QD (circles) and a single MI (stars) with spin 
5 = 1/2 and no nuclear spin. The inset shows the A-shape 
three level optical resonance of bright- and dark-exciton {Xb 
and Xd). The Rabbi-oscillation between Xt, and Xd occurs 
because of the exchange interaction between the exciton and 
the system of MI and nuclear spins. The optical selection rule 
allows decay of X b to vacuum, however, the population of Xd 
decreases indirectly through the conversion of Xd to Xb . 



dark-exciton, is through a transition to bright-exciton 
and spin-flip of MI. This process is schematically de- 
picted in the inset of Fig. [3| A strong dependence of 
dark-exciton population on S is seen in Fig. [2j 

Consistent with the time-evolution of the density ma- 
trix, we propose the exciton wave-function, that fits the 
density matrix via p(t) = \tp(t))(ip(t)\: 

|V(*)> = C bt e- rt ^ 2 \X b A) + C&fVl - e- rt |0,t) 

+C H cos(J cm t/h)e^ 2 \X b , I) + C ol (t)Vl-e- rt \0, |) 

+C dt sm(J cm t/h)\X d ,t)+C di (t)\X d ,i), (4) 

with Cbf = 1/V2. Note that \X b ,i) and \X d ,f) coher- 
ently oscillate because J cm in off-diagonal elements of 
H xm mix these two states. Also from pe6(Tt >> 1) — > 
1/2 we deduce \C ol (Tt » 1)| -> l/%/2, and finally 
C d i(t) = because /O44 = 0. The rest of coefficients 
in \ifi(t)) can be determined numerically by fitting to the 
solutions of density matrix that also fulfill the normal- 
ization of wave- function (tp\tp) = 1. 

A. quantum jump Algorithm 

As noted in Ref. [271 detection of photons requires 
spontaneous emission due to vacuum fluctuations, i.e., 
the photo-emission is a stochastic process, described by 



quantum jump approac H 26 * 27 !. The time evolution of 
the density matrix interacting with photons is given by 
Eq.([3]). At zero temperature in which ng — 0, it is de- 
scribed by a standard Lindblad master equation (ME) 26 

. . (5) 

Comparing Eq.([5| with Eq.([3]) we identify quantum jump 
operator Pm % — |0, M z )(X b , M z \ and its Hermitian con- 
jugate PL = \X b , M z )(0, M z \ that project the excitonic 
state onto vacuum and vice versa without flipping spin of 
MI, hence J2 Mz = \X b ,t)(X b ,t \ + \X b ,i){X b ,i 

I, and {. . .} is the anti-commutator. 

The formulation of quantum jump starts from Eq.([5j). 
In the absence of MI, a recipe for quantum jump al- 
gorithm can be found in Ref. 26 For completeness 
of our presentation first we review this algorithm and 
then generalize it for a system in the presence of MI 
and nuclear spins. We consider optical transition in a 
two-level system consisting of bright-exciton and vac- 
uum without considering an intermediate transition to 
dark-exciton. This condition is fulfilled if we disregard 
presence of any MI and nuclear spin. Here the quan- 
tum jump operators are P — \Q){X b \, P^ = \X b )(0\, 
hence Ptp = \X b )(X b \. Starting at t = with the 
initial condition \ip(t = 0)) = \X b ), we calculate the 
time evolution of the system in a discrete time-steps 
St. In each time-step we evaluate the quantum jump 
probability by calculating 5q Q = T(8t){i/)\P^ P\tp) = T5t 
and drawing a random number r. If r < 5qq a quan- 
tum jump occurs and collapses to |0), otherwise 

\ip(t + St)) = e- rSt( - ptp V 2 \iP{t)). Here the generator for 
the time-evolution operator is a non-Hermitian Hamil- 
tonian H cS = -ihT{P^P)/2. So at t = + St we 
have |V(0 + St)) = e - rSt ^( x ^y 2 \X b ) = er r5t / 2 \X b ) + 
VT— e _r<5 *|0). The last term keeps the norm of \tp) con- 
stant (if we use the norm of wave- function as a constraint 
in our calculation). At this time Sq\ = T(St)e'~ rSt . We 
draw r and if r < Sqo + Sqi then \X b ) — > |0) and a 
photon is detected and calculation is terminated. Oth- 
erwise, \ip{5t + St )) = e -rst (\x b ){x b \)/2^( Q + §t y = 

e~ r( - 2St y 2 \X b ) + Vl-e- r ( 2 *')|0). In nth-step Sq n = 
T(St)e- nrSt thus we calculate an accumulative quantum 
jump probability Sp n = X)fe=o F(5t)e~ krSt , hence 

1 _ e -(n+l)r<5t 

5p n = T(St) 1 _ e _ rgt , (6) 

and if r < Sp n quantum jump occurs. As time 
advances, the chance for a quantum jump becomes 
more likely, however, the probability amplitude for 
X b in \tp) vanishes with the same rate simultane- 
ously. In nth-step if there is still no quantum 
jump, then \ip{nSt )) = e ~ TSt{ -\ x -)^/ 2 \^{[n - l]St)) = 
e- rinSt ^ 2 \X b ) +Vl- e- r (nst) | ) _ 

The quantum jump algorithm in the presence of MI is 
similar to the absence of MI, with a difference that the 
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time-evolution of wave-function is generated by effective 
Hamiltonian H c h = Hqd — ihF6t(P* P)/2 that allows an 
intermediate transition to the dark-exciton due to spin- 
exchange with MI. Therefore the description of quantum 
jump process in the presence of MI is based on a three 
level system depicted in the inset of Fig. |1J and consist 
of |0), \X b ), \X d ) and MI. 



IV. DYNAMICAL EVOLUTION OF MI BY 
TRAIN OF EXCITONS IN THE PRESENCE OF 
NUCLEAR SPINS 

As illustrated in Fig.l a small quantum dot is con- 
tinuoulsy refilled by a non-resonant circularly polarized 
CW laser. The spin polarized excitons transfer into the 
QD containing the complex spin system MI. We assume 
therefore a train of incoming bright excitons \X b ) = | J|, t 
) interacting with MI in the quantum dot. Each electron 
in the exciton transfers spin to MI and creates a superpo- 
sition of dark and bright excitons entangled with nuclear 
spins. At the bright exciton recombination time, t r , pho- 
ton is detected, quantum jump takes place, dark exciton 
wavefunction is erased and spin is transferred from bright 
exciton into a MI and nuclear spin complex and a new 
exciton is created. The exciton removal is performed by 
using the quantum jump projector methocp£12ZI described 
below and yields the wavefunction of the MI and nuclear 
spins. 

The basis for combined exciton-spin sys- 
tem is composed of three groups of basis 
states: vacuum |0, M z , I zl , . . . , I z jy b ), bright ex- 
citon \Xb, M z , I z i, . . . , I Z N b ) and dark exciton 
\X d ,M Zl I z i, . . . ,I zNb ). Only the vacuum and 

bright exciton group of states are coupled to the 
photon field via projectors P A = |0, A)(A b ,A| with 
states | A) = \M z ,I z i, . . . ,I Z N b ) describing a total of 
N s = (2M Z + 1)(2I Z + 1) N " complex spin MI and nuclear 
spin states. In the following we use symbols |A) and \p) 
to represent \M Z , I zX , . . . , I z n„)- 

The time evolution of the density matrix p = \^)(^f\ 
m ME, Eq.((5j), can be generalized by M z — > X. As de- 



scribed in section III A the wave-function 1^) subjected 
to stochastic "birth-death" process of recombination and 
photo-excitation can be used for the time propagation of 
the system coupled with radiation-field and undergoing 
quantum jump process. At t = we start with the initial 
state |*™ =0 (f = 0)) = |0) |xo) of MI and the nuclear spin- 
bath. The index n counts the number of quantum jump 
events. The state = Sa^a = °I^) 1S a ran dom linear 



combination of all possible configurations with the coeffi- 
cients C A being uniformly distributed random complex 
numbers. We note that if we were to compute expecta- 
tion value (M z ) for this random state we would obtain a 

finite value. However, averaging over many sets of 
yields no initial magnetization. 

At t = + a bright exciton created in neighboring QD 
enters the central QD. The creation of |A b ) and annihi- 



lation of |0) are described by operator |A b )(0|, hence the 
wave-function of the system with one exciton is given by 



71=1 



(t = + )) = \X b )(0\V n ="(t = 0)) = \X b )\ Xo ) 



A 

A 



CF 1 (t = 0+)\X b ,\) 



(7) 



where C£ =1 (* = 0+) = C£ =0 (i = 0). The initial wave- 
function of the injected bright exciton, MI and nuclear- 
spins is an uncorrelated state. However, the Hamiltonian 
Hqd that accounts for the exchange coupling, creates 
quantum correlation in the exciton-MI-nuclei complex 
and \^f) evolves into a linear combination of all configu- 
rations, including an entangled state between bright and 
dark excitons. As function of time, the bright-exciton 
decays to vacuum because of coupling with quantized 
electromagnetic-field. 

To be consistent with the quantum jump algorithm 
we discretize the time t into small steps. Note that be- 
cause of the small eh- and Ml-nuclear-spin couplings (J ne , 
J n h and </ n h) the excitons and MI evolve in a frozen- 
fluctuating field of nuclear-spins^"^. The eh recombi- 
nation time, t r , is the smallest time-scale in our model, 
hence St « t r . 

The time evolution of the wave-function of ch-MI- 
nuclei complex is calculated in the Schrodinger picture^ 
by using the relation \^(t+St)) = cx_p(-iH cS 5t/h)\^(t)}. 
Here, is the wave-function of the entire system, 

H c g = Hqd — ih(T/2) ^ A P\P\ where the last term de- 
scribes the decay of bright exciton due to coupling with 
photon-field, hence 

\* n=1 (St)) = cxpi-iHcgdt/h)^ 1 ^ = 0+)). (8) 



In Eq. 



we use exp(—iH c ffSt/h) 



exp(-iH QD 8t/n) exp(-TSt/2 £ A Pi Pa 



0((5t) 2 ) 



with P\P X = \X b , X)(X b , A | and the follow- 
ing identities exp(— T8t ^2 X PjP A )|X b , p) = 
n A exp(-r^P A t P A )|A b , A1 ) = exp(-r<ft)|A%/i), as 
exp(-r<5iP A t F A )|A b ,A) = exp{-T5t)\X b , A), and 
exp(-r^P A t P A )|A b ,/i) = \X b ,n) where fj, ^ A, 
as well as exp(-r<ftPjP A )|X d , A) = \X d , A), and 
exp(-r<5£P;[Px)|0, A) = |0, A). Because H cS is time- 
independent, we employ the method based on Bessel- 
Chebyshev polynomial expansion^ to calculate the time 
evolution of the wave-function 

\y n=1 {5t)) « e-^'E/l^e-s^^'ir^ft = 0+)) 
= e-^^ P ^y\C- x: U6t)\X b ,p) 



C^(5t)\X d ^)}, 



(9) 



Note that Hqd is Hermitian and thus exp(—iHQD5t/h) 
is a unitary operator that conserves the norm of wave- 
function, hence \C^(St)\ 2 + \C^(St)\ 2 = 1. This 
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is in contrast with the operator exp(— TSt/2 -pJ-Px) 
that is non-unitary and does not preserve norm of wave- 
function, however, because it describes the decay of X b 
to vacuum, we build a norm-conserving wave-function by 
adding vacuum. After applying exp(— T8t/2 PlP\) in 
Eq. [9] we find 

l*^ 1 ^*)) = E mM-™/2)C^ x (6t)\X b ) 



-y/l-exp(-TSt)C%f(5t)\0) 
-Cl=\{St)\X d )). 



(10) 



In the limit of V = the coefficients with and with- 
out tilde used in Eqs. [9fl0| are identical. Matching the 
initial conditions between Eq.([7]) and Eq.(10l implies 



CT b \{St = 0) = Cr\t = 0+) and C^ x {St = 0) = 
Using an iterative procedure to propagate the time we 
find 



A 



°0,A 



(*)|0) 



(ii) 



where the coefficients C£=^(t), C^\t), and C x =\(t) are 
determined numerically. This wave-function describes a 
correlated state of bright and dark excitons as well as 
vacuum. Because of the spin flip-flop process of electron 
with MI and nuclear-spins the initially formed bright ex- 
citon \X b ) mixes with the dark-exciton \X d ). 

From the Lindblad ME, the quantum jump transition 
rate is given by rj ump = Tp b with p b representing the 
population of the bright exciton obtained from the full 
QD density matrix after tracing over MI and nuclear 
spin degrees of freedom. The quantum jump probability 
<5f»jump = rj ump <5i is then calculated and compared with 
a random number r generated between zero and one. If 
r < f* r dt5pj ump a quantum jump takes place, photon 
is recorded and the quantum dot is in the ground state. 
The elapsed time t r recorded for this quantum jump is 
the eh-recombination time. At t = t r we allow exciton 
to annihilate by spontaneous emission of a photon. The 
operator that allows annihilation of bright exciton and 
creation of vacuum is 1 0) (X b | ■ Hence 



I*' 



l (t = t r )) = \0)(X b \* n = 1 (t = t r )) 

= |0)E C '^A(* = *r)|A) 
A 

=j2cr 2 (t=t r )\o,x) (12) 



where C{ l=2 (t = t r ) = C x =\(t = t r ). 

Immediately after annihilation of exciton, a new spin 
polarized exciton tunnels into the quantum dot from the 
neighboring dots. The spin polarized exciton interacts 
with the MI spin M and nuclear spins / in a state mod- 
ified by the previous exciton. At t = t r + + , second 
bright exciton X b created in the neighboring QD tunnels 



into the central QD 

|*™ =3 (; = £ r + o+)) = \x b )(o\v n = 2 (t = t r )) 

= \x b )j2cr 2 (t = t r )\\) 

A 

= J2 c r 3 (t = t r +o+)\x b ,\) (13) 



with matching the initial conditions that requires 
C r x l=3 (t = t r + 0+) = C^= 2 (t = t r ). Note that this 
state is not correlated. The quantum correlation appears 
from the time evolution of wave-function generated by 
exchange couplings in H c g right after t = t r 

\^ 3 (t)) = CXp(-lH cB t/h)\t> n = 3 (t = t r +0 + )) 



Ei A )[^r 3 AWi^>+Q,i 3 wio) 



+CT\(t)\x d )}, 



(14) 



where CVr\{t 



x b M"r - 0+) = C'r 3 (tr + 0+) and C x = 3 (t r 



+ ) = are matching conditions. As we see, there is no a 
type of linear combination between |0) and {|A{,), |A^)}, 
because there is no Rabi-oscillations between vacuum and 
excitons. 

To summarize the above procedure and make connec- 
tion across tunneling of exciton and photo-emission we 



formally introduce a projector Q r 



in nth step of quantum jump. The superscripts re- 
fer to annihilated nth and created n + 1th exciton. 
Note that Q„^„+i|A b ") = |A fc " +1 ) and Q n ^ n+1 \X d ) = 
Qji-hi+i|0) = 0, hence the quantum jump operator 
projects out any correlated state composed of superpo- 
sition of bright and dark exciton to a new born bright 
exciton. The new wavefunction in the QD then can be 
constructed as \^f(t = = Q n ^n+i\^{t — t~)) where 
tf = t r ± rj and r\ — > 0. In this state, \X^ +1 ) is initially 
uncorrelated from MI and nuclear-spins. At t = tj it can 
be expressed as 



= | 



(15) 



We observe that after quantum jump the new injected 
exciton X^ +l starts with the normalized (factor An) 
state of MI and nuclear spins J2\ a(^)I^) which was 
left over by the previous bright exciton X£ at the time 
of radiative recombination. Detecting a photon erased 
the dark exciton wave-function and modified the state of 
both MI and nuclear spins. This is the DMP mechanism 
discussed here. With initial condition established, the 
time evolution of the entangled state of photo-carriers 
with MI and spin-bath then can be calculated after up- 
dating the coefficients C's. At the end one needs to av- 
erage over initial conditions. Although the procedure 
discussed here describe the immediate refilling of central 
QD right after annihilation of the exciton, we can im- 
plement a waiting time between the recombination and 
refilling process. 
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FIG. 4: Time evolution of S z of a train of injected photo- 
electrons inside QD (circles), a single MI (stars) and the av- 
erage of Nt = 15 nuclear-spin polarization (triangles). 
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FIG. 5: Time evolution of ensemble average of 20 trajectories 
for S z of a train of injected photo-electrons inside QD (circles) , 
a single MI (stars) and the average of Nt = 15 nuclear-spin 
polarization (triangles) . 



V. NUMERICAL RESULTS AND DISCUSSION 

Our approach to DMP is illustrated for parameters 
based on (Cd,Mn)Te QDs with J cm = 15 meV nm 3 , 
Jhm = 60 meV nm 3 corresponding to the exchange cou- 
pling in the bulk materials, hence J cm = J O m|0e(-Rm)| 2 
and Jhm = Jhm|0h(-Rm)| 2 - The symmetrical quantum 
dot is assumed with Ai = 0. Here </> c /h(-R m ) is the c/h 
envelop- wavefunction in the central dot at R m , the posi- 
tion of MI. J h = 0.6 me\^2l, and J no , J n h, im, and J nn > 
are initialized as random numbers with a mean value in 
the order of 1 /xeV, however, the realistic value for nu- 
clear hyperfine interaction is reported within 1 neV^Sl, 
three orders of magnitude smaller than the energy scales 
used in our finite size calculation. 

Here we discuss numerical results with nuclear spins, 
immediately after refilling of QD by bright exciton. 
Figs. |4j [5] and |6] illustrate DMP/DNP and the quan- 
tum jump trajectories for exciton, MI and nuclear-spins. 
In Fig. [4] a single quantum jump trajectory for spin-1/2 
MI is plotted. In Figs. [5] and [6] the ensemble average of 
twenty quantum jump trajectories for DMP/DNP by a 
series of the injected photo-carriers for spin-1/2 (Fig. [5]) 
and spin-5/2 Mn (Fig. |6| are plotted. Each curve con- 
sists of thousands of time-steps and points. To assign the 
legends, every hundred points, symbols like circle, star 
and triangle are superimposed on each curve. As shown 
the MI and the average polarization of Nj, = 15 nuclear 
spins gradually builds-up by a train of injected bright ex- 
citons. At t = t r , one pair of eh collapses to vacuum with 
(S e z (t)) < 1/2 as part of the e-spin is transferred to MI. 
An empty dot instantaneously absorbs the second photo- 



generated eh pair with total angular momentum j z = — 1 
that transfers to spin of MI and nuclei before its removal. 
We repeat this procedure until the spin polarization in 
MI and nuclear-spins builds-up. The method presented 
here is limited to a finite number of nuclear spins because 
of exponentially increasing with the number of nuclear 
spins computational effort. However, a systematic study 
of the convergence of the numerical results by increasing 
Nf, shows satisfactory outcomes around Nf, = 15. For 
Mn, each of five Mn electrons interact with the exciton 
and nuclear spins with the same exchange coupling be- 
cause of a symmetry that Mn electrons follow to interact 
with the environment. As it is shown in Fig. [6j after ap- 
proximately ten injection and annihilation of excitons in 
the central QD, (S z ) approaches to 5/2 for Mn, reaching 
the maximum polarization. 



VI. SUMMARY 

In conclusion, dynamical magnetic and/or nuclear po- 
larization in complex spin systems is discussed on the 
example of transfer of spin from exciton to the central 
spin of magnetic impurity in a quantum dot in the pres- 
ence of a finite number of nuclear spins. The exciton is 
described in terms of the electron and heavy hole spins 
interacting via exchange interaction with magnetic impu- 
rity, via hyperfine interaction with a finite number of nu- 
clear spins and via dipolc interaction with photons. The 
time-evolution of the exciton, magnetic impurity and nu- 
clear spins is calculated exactly between quantum jumps 
corresponding to exciton radiative recombination. The 




FIG. 6: Time evolution of ensemble average of 20 trajectories 
for S z of a train of injected photo-electrons inside QD (circles), 
a single spin-5/2 Mn with five localized d-electrons (stars) and 
the average of Nb = 15 nuclear-spin polarization (triangles). 



collapse of the wave- function and the refilling of the quan- 
tum dot with new spin polarized exciton is shown to lead 
to build up of magnetization of the magnetic impurity as 
well as nuclear spin polarization. The competition be- 
tween electron spin transfer to magnetic impurity and to 
nuclear spins simultaneous with the creation of dark ex- 
citons is elucidated. The technique presented here opens 
up the possibility of studying optically induced Dynam- 
ical Magnetic Polarization in Complex Spin Systems. 



VII. ACKNOWLEDGEMENT 



RMA thanks Steve Girvin for useful discussion and 
the support from Texas Advanced Computing Center 
(TACC) for computer resources. 



1 R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, 
and L. M. K. Vandersypen, Rev. Mod. Phys. 79, 1217 
(2007). 

2 M. Ciorga, A. Wensauer, M. Pioro-Ladriere, M. Korkusin- 
ski, J. Kyriakidis, A. S. Sachrajda, and P. Hawrylak, 
Phys. Rev. Lett. 88, 256804 (2002); M. Pioro-Ladriere, 
M. Ciorga, J. Lapointe, P. Zawadzki, M. Korkusinski, 
P. Hawrylak, and A. S. Sachrajda, Phys. Rev. Lett. 91, 
026803 (2003). 

3 J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. 
Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, and 
A. C. Gossard, Science 309, 2180 (2005); C. F. Destefani, 
C. McDonald, R. M. Abolfath, P. Hawrylak, T. Brabec, 
Phys. Rev. B 78, 165331 (2008). 

4 M. Bayer, P. Hawrylak, K. Hinzer, S. Fafard, M. Korkusin- 
ski, Z. R. Wasilewski, O. Stern and A. Forchel, Science 
291, 451 (2001). 

5 S. T. Ochsenbein, Y. Feng, K. M. Whitaker, E. Badaeva, 
W. K. Liu, X. Li, and D. R. Gamelin, Nature Nanotech- 
nology 4, 681 (2009); R. Beaulac, Y. Feng, J. W. May, 
E. Badaeva, D. R. Gamelin, and X. Li, Phys. Rev. B 84, 
195324 (2011); K. M. Whitaker, M. Raskin, G. Kiliani, K. 
Beha, S. T. Ochsenbein, N. Janssen, M. Fonin, U. Riidiger, 
A. Leitenstorfer, D. R. Gamelin, R. Bratschitsch, Nano 
Letters, 11, 3355 (2011). 

6 H. O. H. Churchill et al., Nature Phys. 5, 321 (2009); H. 
O. H. Churchill, F. Kuemmeth, J. W. Harlow, A. J. Best- 
wick, E. I. Rashba, K. Flensberg, C. H. Stwertka, T. Tay- 
chatanapat, S. K. Watson, and C. M. Marcus, Phys. Rev. 
Lett. 102, 166802 (2009); S. J. Chorley, J. Wabnig, Z. V. 
Penfold-Fitch, K. D. Petersson, J. Frake, C. G. Smith, and 
M. R. Buitelaar Phys. Rev. Lett. 108, 036802 (2012). 

7 L. C. Bassett, F. J. Heremans, C. G. Yale, B. B. Buck- 



ley, and D. D. Awschalom, Phys. Rev. Lett. 107, 266403 
(2011); G. D. Fuchs, G. Burkard, P. V. Klimov and D. D. 
Awschalom, Nature Physics 7, 789 (2011). 

8 M. Goryca et al. Phys. Rev. Lett. 103, 087401 (2009); 
Phys. Rev. Lett. 102, 046408 (2009). 

9 C. Le Gall et al, Phys. Rev. Lett. 102, 127402 (2009). 

10 L. Besomebes et al., Solid State Commun. 149, 1472 
(2009); L. Besombes, Y. Leger, L. Maingault, D. Ferrand, 
H. Mariette, J. Cibert, Phys. Rev. Lett. 93, 207403 (2004); 
Y. Leger, L. Besombes, J. Fernandez-Rossier, L. Main- 
gault, H. Mariette, Phys. Rev. Lett. 97, 107401 (2006); 
Y. Leger, L. Besombes, L. Maingault, D. Ferrand, H. Ma- 
riette, Phys. Rev. Lett. 95, 047403 (2005); L. Besombes, Y. 
Leger, L. Maingault, D. Ferrand, and H. Mariette, Phys. 
Rev. B 71, 161307 (2005). 

11 S. Mackowski, T. Gurung, T. A. Nguyen, H. E. Jackson, 
L. M. Smith, G. Karczewski, and J. Kossut, Appl. Phys. 
Lett. 84, 3337 (2004). 

12 C. Gould, A. Slobodskyy, D. Supp, T. Slobodskyy, P. 
Grabs, P. Hawrylak, F. Qu, G. Schmidt, and L. W. 
Molenkamp, Phys. Rev. Lett. 97, 017202 (2006). 

13 L. Klopotowski, L. Cywinski, P. Wojnar, V. Voliotis, K. 
Fronc, T. Kazimierczuk, A. Golnik, M. Ravaro, R. Grous- 
son, G. Karczewski, and T. Wojtowicz, Phys. Rev. B 83, 
081306 (2011). 

14 R. Viswanatha, J. M. Pietryga, V. I. Klimov, and S. A. 
Crooker, Phys. Rev. Lett. 107, 067402 (2011). 

15 E. Baudin, E. Benjamin, A. Lemaitre, and O. Krebs, Phys. 
Rev. Lett. 107, 197402 (2011). 

16 J. K. Furdyna, J. Appl. Phys. 64, R29 (1988). 

17 J. Fernandez-Rossier and L. Brey, Phys. Rev. Lett. 93, 
117201 (2004). 

18 A. O. Govorov, Phys. Rev. B 72, 075359 (2005). 



9 



F. Qu and P. Hawrylak, Phys. Rev. Lett. 95, 217206 
(2005); ibid 96, 157201 (2006). 

R. M. Abolfath, A. Petukhov, I. Zutic, Phys. Rev. Lett. 

101, 207202 (2008); R. M. Abolfath, P. Hawrylak, I. Zu- 
tic, Phys. Rev. Lett. 98, 207203 (2007); New Journal of 
Physics 9, 353 (2007); R. M. Abolfath, Phys. Rev. B 80, 
165332 (2009). 

S.-J. Cheng and P. Hawrylak, Eur. Phys. Lett. 81, 37005 
(2008). 

A. Kirilyuk, A. V. Kimel, and T. Rasing, Rev. Mod. Phys. 
82, 2731 (2010). 

A. Oiwa, Y. Mitsumori, R. Moriya, T. Slupinski, and 
H. Munekata, Phys. Rev. Lett. 88, 137202 (2002); H. 
Ohno, D. Chiba, F. Matsukura, T. Omiya, E. Abe, T. 
Dietl, Y. Ohno, K. Ohtani, Nature 408, 944 (2000); D. 
D. Awschalom, J. Warnock, and S. von Molnr, Phys. Rev. 
Lett. 58, 812 (1987) 

D. E. Reiter, T. Kuhn, and V. M. Axt, Phys. Rev. Lett. 

102, 177403 (2009). 

L. Cywinski, Phys. Rev. B 82, 075321 (2010). 
H. J. Carmichael, Phys. Rev. Lett. 70, 2273 (1993); M. B. 
Plenio and P. L. Knight, Rev. Mod. Phys. 70, 101 (1998); 
H. J. Carmichael, An Open Systems Approach to Quantum 
Optics, Lecture Notes in Physics (Springer, Berlin, 1993). 
U. Hohenester, S. S. Communic. 118, 151 (2001); U. Ho- 



henester, C. Sifel, P. Koskinen, Phys. Rev. B 68, 245304 

(2003) . 

M. O. Scully, M. S. Zubairy, Quantum Optics, (Cambridge, 
2001). 

M. Bayer et al, Phys. Rev. B 65, 195315 (2002). 
Optical Orientation, edited by B. Meier and B. P. Za- 
kharchenia (North Holland, Amsterdam, 1984), Ch. III. 
I. A. Merkulov, Al. L. Efros, M. Rosen, Phys. Rev. B 65 
205309 (2002); W. Zhang at al, J. Phys. Cond. Matt. 19, 
083202 (2007). 

A. Melikidze, V. V. Dobrovitski, H. A. De Raedt, M. I. 
Katsnelson, and B. N. Harmon, Phys. Rev. B 70, 014435 

(2004) ; V. V. Dobrovitski, H. A. De Raedt, Phys. Rev. E 
67, 056702 (2003). 

R. M. Abolfath and T. Brabec, Phys. Rev. B 82, 075311 
(2010). 

W. A. Coish, and J. Baugh, Phys. Status Solidi B, 246, 
2203 (2009); W. A. Coish and D. Loss, Phys. Rev. B 72, 
125337 (2005); C. Deng and X. Hu, Phys. Rev. B 73, 
241303(R) (2006); ibid. 74, 129902 (2006); W. Yao, R.- 

B. Liu, and L. J. Sham, Phys. Rev. B 74, 195301 (2006); 
L. Cywinski, W.M. Witzel, S. DasSarma, Phys. Rev. Lett. 
102, 057601 (2009); ibid. Phys. Rev. B 79, 245314 (2009). 



